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Abstract 

A new parameterization of the Jacobi inversion problem is used along with the dynamics 
of the peaks to describe finite time interaction of peakon weak solutions of the Shallow Water 
equation. 

1 Introduction 

Camassa and Holm [1] described classes of n-soliton peaked weak solutions, or "peakons," for an 
integrable (SW) equation 

U t + 3UU X = U xxt + 2U X U XX + UU XXX - 2kU x , (1.1) 

arising in the context of shallow water theory. Of particular interest is their description of peakon 

dynamics in terms of a system of completely integrable Hamiltonian equations for the locations of 
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the "peaks" of the solution, the points at which its spatial derivative changes sign. (Peakons have 
discontinuities in the ^-derivative but both one-sided derivatives exist and differ only by a sign. 
This makes peakons different from cuspons considered earlier in the literature.) In other words, 
each peakon solution can be associated with a mechanical system of moving particles. Calogero [2] 
and Calogero and Francoise [3] further extended the class of mechanical systems of this type. 

For the KdV equation, the spectral parameter A appears linearly in the potential of the corre- 
sponding Schrodinger equation: V = u — A in the context of the inverse scattering transform (1ST) 
method (see Ablowitz and Segur [4]). In contrast, the equation (|1.1|), as well as A^-component 
systems in general, were shown to be connected to the energy dependent Schrodinger operators 
with potentials with poles in the spectral parameter. 

Alber et al. [5,6] showed that the presence of a pole in the potential is essential in a special 
limiting procedure that allows for the formation of "billiard solutions" . By using algebraic-geometric 
methods, one finds that these billiard solutions are related to finite dimensional integrable dynamical 
systems with reflections. This provides a short-cut to the study of quasi-periodic and solitonic 
billiard solutions of nonlinear PDE's. This method can be used for a number of equations including 
the shallow water equation ( |l.l|) , the Dym type equation, as well as A^-component systems with 
poles and the equations in their hierarchies [7]. More information on algebraic-geometric methods 
for integrable systems can be found in [8] and on billiards in [9,10,11]. 

In this paper we consider singular limits of quasi-periodic solutions when the spectral curve 
becomes singular and its arithmetic genus drops to zero. The solutions are then expressed in 
terms of purely exponential r-functions and they describe the finite time interaction of 2 solitary 



new parameterization. First a profile of the 2-peakon solution is described by considering different 
parameterizations for the associated Jacobi inversion problem on three subintervals of the X- 




Namely, we invert the equations obtained by using a 
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axis and by gluing these pieces of the profile together. The dynamics of such solutions is then 
described by combining these profiles with the dynamics of the peaks of the solution in the form 
developed earlier in Alber et al. [9,10]. This concludes a derivation in the context of the algebraic 
geometric approach of the n-peakon ansatz which was used in the initial papers [1,12] for obtaining 
Hamiltonian systems for peaks. More recently n-peakon waves were studied in [13] and [14]. 

The problem of describing complex traveling wave and quasi-periodic solutions of the equation 



(1.1) can be reduced to solving finite-dimensional Hamiltonian systems on symmetric products of 
hyperelliptic curves. Namely, according to Alber et al [5,6,7], such solutions can be represented in 
the case of two-phase quasi-periodic solutions in the following form 

U(x,t) =m + H2~M, (1.2) 

where M is a constant and the evolution of the variables \i\ and \x<i is given by the equations 

dt k = l. 



(1.3) 

dx k = 2. 



Here R([a) is a polynomial of degree 6 of the form R{(x) = M 11^=1 mi). The constant from ( |1.2[) 
takes the form M = 1/2 ^ . Notice that (|1.3| ) describes quasi-periodic motion on tori of genus 
2. In the limit m\ — > 0, the solution develops peaks. (For details see Alber and Fedorov [7].) 

Interaction of Two Peakons. In the limit when ni2 — > m% — > a\ and —* m§ — * 02, we have 
2 solitary peakons interacting with each other. For this 2 peakon case, we derive the general form 
of a profile for a fixed t (t = to, dt = 0) and then see how this profile changes with time knowing 
how the peaks evolve. Notice that the limit depends on the choice of the branches of the square 
roots present in (|1.3| ) meaning choosing a particular sign lj in front of each root. The problem of 
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finding the profile, after applying the above limits to Ql.3|) gives 



h — -, z+h — -, r = a 2 =a 2 dY (1.4) 

/xi(/ii-ai) M2(M2-ai) 

h — 7 r+t2 — -, r = ai = aidr (1.5) 

/xi(^i-a 2 ) /U 2 (/x2-a 2 ) ^1^2 



where 1" is a new variable. This is a new 'parameterization of the Jacobi inversion problem (1.3) 
which makes the existence of three different branches of the solution obvious. In general, we 
consider three different cases: (l\ = 1, l 2 = 1), (h = 1, Z2 = — 1) an d (h = — 1, Z2 = ~~ !)• I n each 
case we integrate and invert the integrals to calculate the symmetric polynomial + fi 2 ). After 
substituting these expressions into the trace formula ( |1.2| ) for the solution, this results in three 
different parts of the profile defined on different subintervals on the real line. The union of these 
subintervals gives the whole line. On the last step these three parts are glued together to obtain a 
wave profile with two peaks. 

The new parameterization dX = ni^dY plays an important role in our approach. In what 
follows each Hi(Y) will be defined on the whole real Y line. However, the transformation from Y 
back to X is not surjective so that ^i(X) is only defined on a segment of the real axis. This is why 
different branches are needed to construct a solution on the entire real X line. 

In the case (h = l 2 = 1), if we assume that there is always one \i variable between a\ and a 2 
and one between and a\ and that initial conditions are chosen so that < < ai < fJ- 2 < a 2 , 
then we find that: /Ui + fi 2 = a\ + a 2 — (mi + ni)aia 2 e x . This solution is valid on the domain 

X < - log(aini + a 2 mi) = X{ , 

where n\,m\ are constants depending on Hi,n 2 . At the point X^ , 

a 2 m\ + a\ni 

Now we consider {l\ = —l,l 2 = 1). Here we find the following expression for the symmetric 



polynomial 

(a 2 - a 1 )e~ x + m 2 n 2 (a 2 - ai)e x 

Hi + Hi = a t + a 2 ■ , 

m 2 + n 2 

which is only defined on the interval 

n 2 a 1 + m 2 a 2 ^ v ^ . a 2 - a x „ + 

log r > X > log ■ = X{ . 

m 2 n 2 (a 2 — a\) m 2 a\ + n 2 a 2 

m 2 ,n 2 are constants which must be chosen so that both hi and \i 2 are continuous at and that 

the ends of the branches match up, that is so that X^ = X± . These conditions are satisfied if 

m 2 = —(a 2 -ai)mi, (1.6) 
ai 

n 2 = —{a 2 -ai)ni. (1.7) 
Continuing in this fashion we arrive at the final 3 branched profile for a fixed i, 

U = -{aiM + a 2 N)e x if X < - \og{N + M) (1.8) 

ai a 2 e- x + M N e x {a 2 -a 1 ) 2 

a 2 M + aiN ^ ' ' 

-W + M)<X<* ^ + 4» N (1.10, 

where we have made the substitution M = a 2 m\ and N = a\ri\ and used the trace formula (|1.2|), 
Please place the first figure near here. 



Time evolution. So far only a profile has been derived. Now we will include the time evolution 
of the peaks to find the general solution for the two peakon case. To do this we use functions qi(t) 
for i = 1, 2 introduced in Alber et al. [9] 

m(x = qi(t),t) = o, 

for all t and i = 1, 2 which describe the evolution of the peaks. All peaks belong to a zero level set: 
Hi = 0. Here the /i-coordinates, generalized elliptic coordinates, are used to describe the positions 



of the peaks. This yields a connection between x and t along trajectories of the peaks resulting in 
a system of equations for the The solutions of this system are given by 

qi (t) = ql-a 2 t-\og\l-C l e^- a ^ t \+\og{l-C l ) (1.12) 
q 2 (t) = q 0-a 2 t + log\l-C 2 e^- a ^\-log(l-C 2 ), (1.13) 

where Q = (^(0) - ai)/(^(0) - a 2 ). 

The solution defined in ( |l.8[) has the peaks given in terms of the parameters N and M. So to 
obtain the solution in terms of both x and t, these parameters must be considered as functions of 
time. The complete solution now has the form 

U = -(oiM(t) +a 2 N(t))e x if X < - log(iV(t) + M{t)) (1.14) 
aia 2 e~ x + M(i) iV(i) e x {a 2 - m) 2 



U 



a 2 M(t) + ai N(t) 

4M(t) + a 2 N(t) 



if - log (tf(t) + M(t)) < X < log ^r^^ (I-") 

x q|M(t) + Q fJV(t) alM(t)+ajJV(t) 
M(t) JV(t)(o2 - ai ) 2 ^ S (aa - 0l ) 2 M(t) JV(t) ' 1 j 

where the functions M(t),N(t) are determined by the relations 

AT(t) + M(t) = e -*W = - 11 ^° 16 1 (1.17) 

1 — Ci 

a|M(t) + afiV(t) _ 2 ffl _ (a 2 - Ql ) - C 2 e~^\ 

M{t) N{t) ~ 12 lj (1 - C 2 ) ' 1 j 

where <zi(i), 92 (^) are taken from ( |1.12| )- (|1.13| ), This system can be solved to find that 

a 2 - a\ + ^(t)g(t) ± ./(a 2 - a 2 ) 2 - 2A(t)B{t)(a{ + a 2 ) + A(t) 2 B(t) 2 
[t) ~ 2B(t) 1 j 

N(t) = A(t)-M(t), (1.20) 
where A(t) = e~ qi ^ and -B(i) = (a 2 — a\) 2 e q2< ^ > . These functions contain 4 parameters, but in fact 
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these can be reduced to two parameters by using the following relations 

, l( o) . - log <M ( o> rf ( ,)- sg||±agp ( , 21) 

, n , ,„„ a|M(0) + afiV(O) a^^MfO) + fll iV(0)) 

92(0) = bg (a 2 - ai ) 2 M (0) W) q2i )= alM(0)+a?JV(0) ' (L22) 

Some care must be used in choosing the sign in (|1.1S| ). It is clear that for large negative t, fxi(qi(t),t) 

refers to the path of one peakon while for large positive t it refers to the other. If this were not the 

case, simple asymptotic analysis of ( |1.12 ) would show that the peakons change speed which is not 



the case. Therefore qi(t) represents the path of one of the peakons until some time t* and the other 
one after this time. The opposite is true for q 2 (t). At the time t* we say that a change of identity 
has taken place, t* can be found explicitly by using the fact that at this time, the two peaks must 
have the same height. But the peaks have the same height exactly when 

a 2 M(t*) = ai N(t*). (1.23) 

Without loss of generality we can rescale time such that t* = 0. In this case ( |1.23| ), due to the 
original definitions of mi,ni given in terms of li® fi 2 , corresponds to a restriction on the choice of 
/J( and fi®, namely 

2/jO-ai 2 ii\- a 2 

~ a 2— = a i— • ( L24 ) 

H\ -a 2 Hi - ai 

This condition is satisfied for example when ll\ = aia2 an d li^ = 1 ° 2 . Also notice that under 

a\ + a 2 2 

this rescaling, the phase shift is simply gi(0) — q 2 (0). 
Please place the second figure near here 

So we now have a procedure to make the change of identity occur at t = 0, i.e. Hi goes from 
representing the first peakon to the second one at t = 0. This change is represented by the change 
in the sign of the plus/minus in ( 1.1S| ). That is, the sign is chosen as positive for t < and negative 



for t > 0. However, M remains continuous despite this sign change since the change of identity 



occurs precisely when the term under the square root is zero. Therefore (p^)-(p^) and (p^| ) 
together describe the solution U(X, t) of the SW equation as a function of x and t depending on 
two parameters M(0), N(0). 

By using the approach of this paper weak billiard solutions can be obtained for the whole class 
of ra-peakon solutions of iV-component systems. 
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Figure 1.1: This is a plot U(x, 0), a profile of the solution to the SW equation for the two peakon 
case where a\ = — l,ct2 = —3,^5 = —.h,^ = —1.3. Notice how the solution is defined on three 
different branches. 
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Figure 1.2: This is a plot of (fii+122), which is what we are seeking in this section. The parameters 
used are a\ = —1, 02 = —3, rf[ = —.5, [i® = —1.3. Again notice how the solution is defined on three 
different branches. 
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